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Abstract. A self-avoiding walk (SAW) on the square lattice is prudent if it never takes a 
step towards a vertex it has already visited. Prudent walks differ from most classes of SAW 
that have been counted so far in that they can wind around their starting point. 

Their enumeration was first addressed by Prea in 1997. He defined 4 classes of prudent 
walks, of increasing generality, and wrote a system of recurrence relations for each of them. 
However, these relations involve more and more parameters as the generality of the class 
increases. 

The first class actually consists of partially directed walks, and its generating function is 
well-known to be rational. The second class was proved to have an algebraic (quadratic) 
generating function by Duchi (2005). Here, we solve exactly the third class, which turns out 
to be much more complex: its generating function is not algebraic, nor even D-finite. 

The fourth class — general prudent walks — is the only isotropic one, and still defeats us. 
However, we design an isotropic family of prudent walks on the triangular lattice, which we 
count exactly. Again, the generating function is proved to be non-D-finite. 

We also study the asymptotic properties of these classes of walks, with the (somewhat 
disappointing) conclusion that their endpoint moves away from the origin at a positive speed. 
This is confirmed visually by the random generation procedures we have designed. 



1. Introduction 

1.1. Families of self- avoiding walks 

The study of self-avoiding walks is a famous "elementary" problem in combinatorics, which is 
also of interest in probability theory and in statistical physics [21]. Recall that, given a lattice 
with some origin O, a self-avoiding walk (SAW) is a lattice path starting from O that does not 
visit the same vertex twice (Fig. [T|). 




Figure 1. A self-avoiding walk on the square lattice, and a (quasi-)random 
SAW of length 1,000,000, constructed by Kennedy using a pivot algorithm |19j . 

It is strongly believed that, for two-dimensional lattices, the number c{n) of n-step SAW and 
the average end-to-end distance of these walks satisfy 

c{n) ^ af/^rP and E(£)„) ~ ku'^ 
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where 7 ~ 11/32 and ly = 3/4. The growth constant fj. is lattice-dependent. Several independent 
but so far non-rigorous methods predict that = \/2 + y/2 on the honeycomb lattice. Moreover, 
numerical studies suggest that ^ may also be a bi-quadratic number for the square lattice |18) . 
On the probability side, it has been proved that, if the scaling limit of SAW exists and has 
some conformal invariance property, it must be described by the process SLE(8/3) (stochastic 
Loewner evolution) [20]. This would imply that the predicted values of 7 and are correct. 

The fact that all these conjectures only deal with asymptotic properties of SAW tells us how 
far the problem is from the reach of exact enumeration. The followers of this discipline thus 
focus on the study of sub-classes of SAW. A simple family consists of partially directed walks, 
that is, self-avoiding walks formed of North, East and West steps. It is easy to see that their 
generating function is rational [2H1 Example 4.1.2], 



l-2t-t^' 



which gives c(n) ^ Q^m") with ^ = 1 + V2 = 2.41... The above series can be refined by taking 
into account the coordinates (X„,y„) of the endpoint, and the analysis of the result gives: 

E(X„) = 0, E(A:2)~n and E(r„) ~ 

We use the notation shorthand for a„ ^ a5„ for some positive constant a. 

The prudent walks studied in this paper form a more general class of SAW which have a 
natural kinetic description: a walk is prudent if it never takes a step pointing towards a vertex it 
has already visited. In other words, the walk is so cautious that it only takes steps in directions 
where the road is perfectly clear. Various examples are shown on Fig. [2l In particular, partially 
directed walks are prudent. Note that it is equivalent to require that the new step avoids the 
past convex hull of the walk: hence prudent walks are the discrete counterpart of the walks 
evolving in discrete time but continuous space studied by Angel et al. [2j. 
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Figure 2. Four prudent walks: the first one is partially directed (or: 1-sided), 
the other ones are respectively 2-sided, 3-sided and 4-sided. 



These walks have already been studied in the past under different names: self-directed 
walks [30], outwardly directed SAW [27], exterior SAW [25], and finally prudent walks [TUfTO]. 
We refrain from the temptation of inventing one more name and stick to the latter terminology. 
The first two papers above deal with Monte-Carlo simulations. Prea [25] was, to our knowledge, 
the first to address enumerative questions. He wrote recurrence relations defining an array of 
numbers c{n;i,j,h) that count prudent walks according to their length (n) and to three addi- 
tional catalytic parameters (i, j, h). By this, we mean that these parameters are essential to the 
existence of these recurrence relations, and that it is far from obvious how to derive from them a 
recursion for, say, some numbers c(n;i, j) that would only take into account two of the catalytic 
parameters (or for c(n)). Prea also defined four natural famihes of prudent walks of increas- 
ing generality, called k-sided, for k ranging from 1 to 4. In particular, 1-sided walks coincide 
with partially directed walks, and 4-sided walks coincide with general prudent walks. (Precise 
definitions will be given below.) He wrote recurrence relations for each of these classes: three 
catalytic parameters are needed for general (4-sided) prudent walks, but two suffice for 3-sided 
walks, while one is enough for 2-sided walks. No catalytic parameter is needed for 1-sided walks. 
This refiects the increasing generality of these four classes of walks. 
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Recall that the generating function of 1-sided walks (partially directed walks) is rational ([T]) . 
Duchi [11] proved that 2-sided walks have an algebraic (quadratic) generating function. She 
also found an algebraic generating function for 3-sided walks, but there was a subtle flaw in her 
derivation, which was detected by Guttmann [16]. He and Dethridge performed a numerical 
study of prudent walks, in order to get an idea of their asymptotic enumeration, and of the 
properties of the associated generating functions [TO]. In particular, they conjectured that the 
length generating function of general prudent walks is not D- finite, that is, does not satisfy any 
linear differential equation with polynomial coefficients. This implies that it is not algebraic. 

1.2. Contents 

In Section [2] of this paper, we collect functional equations that deflne the generating functions 
of the four classes of prudent walks introduced by Prea [25]. This is not really original, as these 
equations are basically equivalent to Prea's recurrence relations. Moreover, similar equations 
were written by Duchi [Tl]. Our presentation may be a bit more systematic. 

In Sections [3] to E] we address the solution of these equations. The case of 1-sided walks is 
immediate and leads the above rational generating function. We then recall how the kernel 
method solves linear equations with one catalytic variable. In particular, it provides the gener- 
ating function of 2-sided walks (Section [31. The extension of this method to linear equations 
with two catalytic variables is not yet completely understood, although a number of papers have 
been devoted to instances of such equations recently [3 O [U [HI [TJl [23l [22] • Underlying the 
bi-variate kernel method is a certain group, which depends on the equation. Roughly speaking, 
the instances that have been studied suggest that the solution is "nice" if the group is flnite. 
This belief is conflrmed by the example of 3-sided walks. The corresponding equation is associ- 
ated with an infinite group, but we can still solve it, and prove that the generating function of 
these walks is not D- finite, having a rather complex singularity structure (Section [1]) . We also 
prove that the growth constants of 3-sided walks and 2-sided walks are the same. It is actually 
predicted from numerical experiments that the growth constant of general prudent walks is also 
the same |10) . 

The final equation, which deals with general prudent walks and involves three catalytic vari- 
ables, still defeats us. This is a bit annoying, as the other classes are by definition anisotropic. 
However, our understanding of the role of catalytic variables leads us to introduce a new isotropic 
class of prudent walks, on the triangular lattice, which are described by (only) two catalytic vari- 
ables (Fig. [31 left). Again, the associated equation corresponds to an infinite group. We solve 
it, and prove that the generating function of triangular prudent walks is not D-finite, having a 
natural boundary (Section [5]). 




Figure 3. Left: A triangular prudent walk in a box of size 7. Right: The two 
catalytic parameters involved in the enumeration. 

We also refine our equations to take into account parameters related to the end-to-end distance 
of a prudent walk: for instance the coordinates of the endpoint, or the size of the smallest 
rectangle containing the walk. Extending our solutions to these refined equations is harmless, 
but the conclusion we draw from these results is somewhat disappointing: the prudent walks we 
can solve drift away from the origin at a positive speed. In other words, the end-to-end distance 
grows linearly with the length of the walk. We do not know what happens for general (4-sided) 
prudent walks. 
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Finally, in Section[6l we address the uniform random generation of n-step prudent walks. Their 
step-by-step recursive structure allows for a standard recursive approach, in the spirit of |24) . 
This approach first requires a precomputation stage, followed by a generation stage, which is 
usually linear. We emphasize that the (costly) precomputation may require less information 
than the generation itself, and we use this to optimize the precomputation stage. Our final 
procedures involve the precomputation and storage of up to O(n^) numbers (for general prudent 
walks), so that the typical length we can reach is a few hundred. This still provides interesting 
pictures (Figs.HlelEllHl). 




Figure 4. Random prudent walks. 

1.3. Families of prudent walks 

Let us conclude this long introduction with some definitions and notations. The box of a 
square lattice walk is the smallest rectangle that contains it. It is not hard to see that the 
endpoint of a prudent walk is always on the border of the box. This means that every new step 
either walks on the border of the box, or moves one of its four edges. 

Using this box, we can give a kinetic description of partially directed walks: a prudent walk 
is partially directed if its endpoint, as the walk grows, always lies on the top of the box. This 
is why partially directed walks are also called 1 -sided. The generalisation of this terminology is 
natural: a prudent walk is 2-sided if its endpoint always lies either on the top, or on the right 
edge of the box. It is 3-sided if its endpoint is always on the top, right or left edge of the bojfl. 
Of course, A-sided walks coincide with general prudent walks (Fig. [H). 

Consider now a walk on the triangular lattice. Define the (triangular) box of the walk as the 
smallest triangle pointing North that contains the walk. The walk is a triangular prudent walk if 
each new step either inflates the box, or walks along one edge of the box in a prudent way (that 
is, not pointing to an already visited vertex). An example is shown in Fig. [3l Note that this is 
not the natural counterpart of a square lattice prudent walk: this counterpart would just require 
the walk to avoid pointing to an already visited vertex. But then the natural "box" would be a 
hexagon: every new step would either inflate the hexagonal box, or walk along its border. As 
we will see below, the number of edge lengths of the box (3 for a hexagon, 1 for a triangle) is 
directly related to the number of catalytic parameters we have to introduce, and this is what 
makes triangular prudent walks relatively easy to handle. 

^As noticed by Uwe Schwerdtfeger, we actually require this to hold as well for non-integral points, that is, 
when the walk is considered as a continuous process. The aim of this is to avoid considering the walk ESW, 
among others, as a 3-sided walk. 
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Given a class of walks C, the generating function of walks of C, counted by their length, is 

wee 

where \w\ denotes the length of the walk w. The generalisation of this definition to the series 
C{t;ui, . . . ,Uk) counting walks according to their length and to k additional parameters is im- 
mediate. We will often drop the length variable t, denoting this series C(wi, . . . ,Uk)- Recall that 
a one-variable series C{t) is algebraic if it satisfies a polynomial equation P[t,C{t)) = 0, and 
D-finite if it satisfies a linear ODE with polynomial coefficients, Pk{t)C^''\t) + ■ ■ ■ + Pi{t)C' (t) + 
Pn{t)C{t) = 0. Every algebraic series is D-finite. We refer to [29] for generalities about these 
classes of power series. 

2. Functional equations 

The construction of functional equations for all the families of prudent walks we study rely 
on the same principle, which we first describe on 1-sided (partially directed) walks. 

Consider a 1-sided walk. If it ends with a horizontal step, we can extend it in two different 
ways: either we repeat the last step, or we change direction and add a North (N) step. Otherwise, 
the walk is either empty or ends with a North step, and we have three ways (N, E and W) to 
extend it. This shows that North steps, which move the top edge of the box, play a special role 
in these walks. Our functional equation is obtained by answering the following question: where 
is the last North step, and what has happened since then? 

More specifically, let P{t) denote the length generating function of 1-sided prudent walks. 
The contribution to P{t) of walks that contain no North step (horizontal walks) is 

^ 1-t 

n>l 

The other walks are obtained by concatenating a 1-sided walk, a North step, and then a horizontal 
walk. Their contribution is thus 

Adding these two contributions gives the equation 

, l + t 1 + t , , 

from which we readily derive the rational expression |(T|). 

The principle of this recursive description extends to /c-sided walks for each k. We say that a 
step of a fc-sided walk is inflating if, at the time it was added to the walk, it shifted one of the k 
edges of the box that are relevant in the definition of fc-sided walks. For instance, when k = 2, 
an infiating step moves the top or right edge of the box. We write our equations by answering 
the question: where is the last infiating step, and what has happened since then? 

Since then, the walk has grown without creating a new inflating step. What does that mean? 
Assume fc > 2, that the last infiating step was North, and that, since then, the walk has taken 
TO East steps. Then to cannot be arbitrarily large, otherwise one or several of these East steps 
would be infiating, having moved the right edge of the box. This observation explains why we 
have to take into account other "catalytic" parameters in our enumeration of prudent walks. 
For instance, for a 2-sided walk, we keep track of the distance between the endpoint and the 
NE corner of the box, using a new variable u (Fig. [5]). For 3-sided walks ending on the top of 
the box, we keep track of the distances between the endpoint and the NE and NW corners of 
the box (variables u and v). For 3-sided walks ending on the right edge of the box, we keep 
track of the distance between the endpoint and the NE corner (variable u) and of the width 
of the box (variable w) . For 4-sided walks ending on the top of the box, we keep track of the 
distances between the endpoint and the NE and NW corners (variables u and v), and of the 
height of the box (variable w). These parameters, and the names of the corresponding variables. 



6 



MIREILLE BOUSQUET-MELOU 



are schematized in Fig. [H They give rise to series with one, two or three catalytic variables. For 
instance, for 4-sided walks ending on the top of their box, we will use the series 



T{t; u, V, w) = T(w, u, w) 



where Tijji = Tij^h{t) counts 4-sided walks ending on the top of their box, at a distance i (resp. 
j) from the NE (resp. NW) corner, such that the height of the box is h. Similar notation will 
be used for the other classes of walks. 



V u 



V u 






J} 




Two-sided 



Three-sided 



Four-sided 



Figure 5. Catalytic variables for /c-sided walks, k >2. 



Finally, for triangular prudent walks ending on the right edge of their (triangular) box, we 
keep track of the distances between the endpoint and the SE and N corners of the box (variables 
u and V, see Fig. [3]). 

2.1. Two-sided prudent walks 

Lemma 1. The generating function T{t]u) = T(u) of 2- sided walks ending on the top of their 
box satisfies 

tu(l-t^) \ , , 1 u~2t , ^ 

1- \ w \A T{u) = ——+t -T{t). 

(l — tu){u — t)J 1 — tu u — t 

The generating function of 2-sided walks, counted by their length and the distance of the endpoint 

to the NE corner of the box, is 

P{t;u) = 2T{t;u)-T{t;0). 

Proof. We partition the set of 2-sided walks ending on the top of their box into 3 classes, 
depending on the existence and direction of the last inflating step (LIS). This step, if it exists, 
has moved the right or top edge of the box. 

1. Neither the top nor the right edge has ever moved: the walk is a sequence of West steps. 
The generating function for this class is 

1 

1-tu' 

2. The LIS goes East. This implies that the endpoint of the walk was on the right edge of 
the box before that step. After that East step, the walk has made a sequence of N steps 
to reach the top of the box. Observe that, by symmetry, the series T{t; u) also counts 
walks ending on the right edge of the box by the length and the distance between the 
endpoint and the NE corner. These two observations give the generating function for 
this class as 

t^T.f = tT{t). 

i>0 

3. The LIS goes North. After this step, there is either an (unbounded) sequence of West 
steps, or a bounded sequence of East steps. This gives the generating function for this 
class as 



-T{u) + t^T,^ t'u^-' = ^^T{^) + Wu) tT{t)) , 



1 — tu ^-^ ^-^ 1 — tu u — t 

i>0 fc=0 
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Adding the 3 terms gives the functional equation satisfied by T{u). 

The expression of P(i; u) rehes on an inclusion-exclusion argument: we first double the contri- 
bution of T(u) to take into account walks ending on the right edge of the box, and then subtract 
the series T(0) counting those that end at the NE corner. B 

2.2. Three-sided prudent walks 

Lemma 2. The generating functions T{t;u,v) = T{u,v) and R(t;u,w) = R{u,w) that count 
respectively 3-sided walks ending on the top and on the right edge of their box are related by 

tuv( 1 — t'^) \ t^v t^u 

1 - ) = 1 + + ^) - - 

tUw(l—t^) \ , ^ 1 ^ t^W , ^ 

1 - Riu, w) = + tT(tw, w) R(t, w). 3 

{u-t){l-tu)J ^ ' ' 1-tu ^ ' ' u-t ' ^' 

The generating function of 3-sided walks, counted by the length and by the width of the box, is 

t 



P{t; u) = T{t; u, u) + 2R(t; 1, u) - 2T{t] u, 0) - 



1 - t 

Proof. We partition the set of 3-sided walks ending on the top of their box into 4 subsets, 
depending on the existence and direction of the LIS, which has moved the right, left or top edge 
of the box. 

1. There is no inflating step at all: the walk is empty and contributes 1 to the generating 
function. 

2. The LIS goes East. This case is analogous to Case (2) of 2-sided walks, with generating 
function 

tv^R,^jfv^ = tvR{t, v). 

id 

3. Symmetrically, the case where the LIS goes West is counted by 

tuR(t, u). 

4. If the LIS is a North step, it is followed by a bounded number of West steps, or by a 
bounded number of East steps. This case is counted by: 



i,j>0 \k=0 k=0 



(uT{u, v) — tvT{tv, v)) H (vT{u, v) — tuT(u, tu)) — tT{u, v). 



u — tv V ~ tu 

Adding the 4 terms gives the first equation of the lemma. 

For 3-sided walks ending on the right edge of their box, the last infiating step cannot go West. 
Three cases remain: 

1. There is no infiating step at all: the walk consists of South steps. The generating function 
for this class is 

1 



l-tu 

2. The LIS goes East. This case is analogous to Case (3) of 2-sided walks: the LIS is 
followed by an unbounded number of South steps, or by a bounded number of North 
steps. The generating function is 

t^uw tw 

-R{u, w) H {uR{u, w) - tR(t, w)) . 



1 — tu u — t 

3. The LIS goes North. This case is analogous to Case (2) of 2-sided walks. The generating 
function is found to be 

tT(tw, w). 
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Adding the 3 terms gives the functional equation for R{u, w). 

The expression of P(t; u) again rehes on an inclusion-exclusion argument, based on the enu- 
meration of walks ending on a prescribed set of edges of their box. ^ 

2.3. General prudent walks on the square lattice 

Lemma 3. The generating function T{t; u, v, w) = T{u, v, w) of prudent walks ending on the 
top of their box satisfies 

(l - ^""^(-'^ — — TT ) T{u,v,w) = 1 + g{w,u) +g{w,v) ^-^^g{v,w) ^-^g{u,w) 

\ [u — tv)[v — tu) J u — tv V — tu 

with Q(u^v) = g(t;u,v) = tvT(t; u,tu,v). 

The generating function of prudent walks, counted by the length and the half-perimeter of the 
box, is 

P{t; u) = 1 + 4T{t; u, u, u) - 4T{t; 0, u, u). 

We have learnt from [lO] that Andrew Rechnitzer has independently obtained the functional 
equation satisfied by T(u, v, w). 

Proof. Again, we partition the set of prudent walk ending on the top of their box into 4 subsets, 
depending on the existence and direction of the last infiating step. Note that the LIS cannot be 
a South step. 

1. There is no inflating step: the walk is empty, and contributes 1 to the generating function. 

2. The LIS goes East. This case is analogous to Case (2) of 2-sided walks. Using the 
obvious symmetry between prudent walks ending on the top and on the right edge of 
their box, we find that the generating function for this class is 

tv ^ T^^jMt^w'+^v'' = tvTiw, tw, v). 

3. Symmetrically, the generating function of prudent walks in which the LIS goes West is 

tuT{w, tw, u). 

4. Finally, the generating function of prudent walks in which the LIS goes North is analo- 
gous to Case (4) of 3-sided walks ending on the top, with generating function: 

wt wt 

{uT(u,v,w) — tvT(tv,v,w)) + (vT(u,v,w) — tuT(u,tu,w)) — twT(u,v,w). 

u — tv V — tu 

Adding the 4 terms provides the functional equation for T{u,v,w), given the obvious symme- 
try T{u,v,w) = T{v,u,w). The expression for P{t;u) again relies on an inclusion-exclusion 
argument, based on the enumeration of walks ending on a prescribed set of edges of their box. ^ 

2.4. Triangular prudent walks 

Lemma 4. The generating function R{t;u,v) = R{u,v) of triangular prudent walks ending on 
the right edge of their box satisfies 

1 - ^ ^/ ^ R(u,v) = l + tu(l+t) R(u,tu) + tv(l + t) R(tv,v). 4 

\ [u — tv)(v — tu) J V — tu u — tv 

The generating function of triangular prudent walks, counted by the length and the size of the 
box, is 

P{t; u) = l + 3R{t; u, u) - 3R{t; u, 0). (5) 

Proof. We partition the set of prudent walks ending on the right edge of their box into 7 subsets, 
depending on the existence and direction of the last infiating step. 

1. There is no infiating step: the walk is empty, and contributes 1 to the generating function. 
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2. The LIS is a NW step. This impHes that the endpoint of the walk was on the left edge 
of the box before that step. Thanks to the obvious symmetry between walks ending on 
the left edge and right edge of their box, we obtain the generating function of walks of 
this type as 

tu ^ it*+-' tuR{tu, u). 

3. The LIS is a W step. Again, the endpoint of the walk was on the left edge of the box 
before that step. The generating function of walks of this type is 

iu^i?,jf+^it'+^' =t^uR{tu,u). 

hi 

The cases where the LIS goes SE or SW are very similar to the two previous cases. The endpoint 
of the walk was on the bottom edge of the box before the last inflating step. 

4. The generating function of walks such that the LIS is a SE step is 

= tvR{v,tv). 

hi 

5. The generating function of walks such that the LIS is a SW step is 

tv J2 R^,jt^^^v'^^ = t^vR{v, tv). 

We are left with the two richer cases where the LIS goes East or North-East. The endpoint of 
the walk was already on the right edge of the box before the last inflating step, and the LIS is 
followed by a bounded sequence of SE or NW steps. 

6. The generating function of walks such that the LIS is an E step is 

^^^1 t''u'-^v'+^ + t'li'+'i'^-M 

i,j \fc=0 fc=l / 

tv t^ u 

= (uR(u, v) — tvRitv, v)) H {vR(u, v) — tuRiu, tu)) . 

u — tv V — tu 

7. The generating function of walks such that the LIS goes NE is 
tu'^Ri^j ^t'^v^-^u 

i,j \fc=0 

t^V 

(vR(u, v) — tuR(u, tu)) H (uR(u, v) — tvRitv, v)) . 

V ~ tu u — tv 

We add the 7 terms above and note that v) = R{v, u) to obtain the functional equation for 
R{u, v). The expression of P{t; u) again relies on an inclusion-exclusion argument. ■ 

3. Enumeration and asymptotic properties of 2-sided prudent walks 

In this section, we recall how the kernel method works on linear equations with one catalytic 
variable [7l[3l[26], using the example of 2-sided walks. We flrst recover Duchi's algebraic generat- 
ing function, and then reflne our enumeration to keep track of other statistics like the end-to-end 
distance of the walk. This section is a sort of warm-up before solving the more difficult equations 
of Sections H] and [H 

Proposition 5. The generating function P{t; u) of 2-sided walks, counted by their length and 
by the distance between the endpoint and the NE corner of the box, is 

2{l~t^)(l-t)U 
^ '""^ {1 - uU){l - tU){2t - U) 
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where 

u ^ u{t) = 1-^+^'+^^- v/a-^^)(i-2^-t^) ^ 

In particular, the length generating function is 



1 - 2i - 2<2 + 2t3 \ ' ' V 1 - 2i - t2 I • 

Proof. We start from the functional equation of Lemma [H written as 

((1 - tu){u ~t)- tu{l - <2)) T{u) =u-t + t{l- tu){u - 2t)T{t). (6) 

The series U = U{t) given in the proposition is the only power series in t that cancels the kernel of 
this equation, that is, the polynomial ((1 — tu){u — t) ~ tu{\ — i^)). The series T{U) = T{t\ U) 
is well-defined. Replacing u by C/ in the equation cancels the left-hand side, and hence the 
right-hand side, giving 



t{l - tU){2t - U) 

Then ^ gives T(u), and the second equation of Lemma [T] provides P[t] u). 



Figure 6. Random 2-sided walks of length 500, 1354 and 3148. 

Proposition 6 (Asymptotic properties of 2-sided walks). The length generating function 
of 2-sided walks has a unique singularity of minimal modulus, p ~ 0.403, which is a simple pole 
and satisfies 1 — 2p — 2p^ + 2p^ ~ 0. The number of n-step 2-sided walks satisfies 

Pn ~ K/i" with a = — ~ 2.48 and n ~ .„ ^ ^ f — ^— rr — 2.51. 

p (3p -I- 1) (5p - 2) 

The distance between the endpoint and the NE corner of the box in a random n-step 2-sided 
walk follows asymptotically a geometric law of parameter 2p. In particular, the average value of 
this distance is asymptotically constant, equal to — 4.15. 

Let Xn (resp. Yn) denote the abscissa (resp. ordinate) of the endpoint of a random n-step 
2-sided walk. Then the mean and variance of Xn -\- Yn satisfy 

E{Xn + Yn)^mn, V(X„ + y„) ~ n, 

where 

m=^c.0.63 and ^ ^(p+lfp ^ 
3/0+1 (3p+lf(l-p) 
and the variable " converges in law to a standard normal distribution. 

Finally, 

EiXn - Yn) = 0, Y{Xn - Yn) ^ 5^ 71, 

where 

p{p'-2)il+p) ^ 
(p2 + p-i)(3p-l)(l + 3p) ^ • ' 
and the variable "^^'y^" converges in law to a standard normal distribution. 

These asymptotic properties are in good agreement with the random 2-sided walks of Fig. [6l 



PRUDENT SELF-AVOIDING WALKS 



11 



Proof. We start from the expression of P(t; 1) given in Proposition [H The singularities of P{t; 1) 
are found among the roots of the denominator (1 — 2i — 2t^ + 2t^) and those of the discriminant 
{l-f^){\-2t-t^). It is not hard to see that the smallest one (in modulus) is p, corresponding to 
a simple pole of the series. In the expression of P{t] u) in terms of U given in Proposition [U this 
pole is found when U = 2t. This implies ^(1 — -2p - p^) = (1 - 4p^)/2. The asymptotic 
behaviour Pn ^ up," easily follows, with p, = p^^. 

To study the distance to the NE corner, consider P{t;u) (given in Proposition [5]) . For u in 
a neighborhood of 1, P{t;u) still admits p as its unique dominant singularity, which remains a 
simple pole. This gives 

r]Pit;u) ^ l-^^[npit;i), 

with U{p) = 2p, and the result follows using a continuity theorem for probability generating 
functions [H Thm. IX. 1]. 

Then, we enrich our enumeration by taking into account the sum of the coordinates of the 
endpoint, using a new variable z. The functional equations of Lemma [1] become: 

tuz^(l — t'^) \ , z u — 2tz^, 

^-1 rV( T-AT{t,z;u) = -+tz —T{t,z;tz), 

[z — tu)[u ~ tz) J z ~ tu u — tz 

P{t,z;u) = 2T{t,z;u)~T{t,z;0). 

We solve them in the same way that we solved the case z — 1, and obtain: 

PU ^ _ 2z^{l-t'Kl-tz)U _ 
^ ' ^' ~ (z2 - uU){z - tU){2tz - U) 

with 



^ l-tz + e+ t^Z - J[l - t:^)[l + t - tz + t2z)(l -t-tz- t^z) 

U = U{t,z) = z ^ . 

For z in a neighborhood of 1, the radius of convergence pz of P{t, z; 1) is reached when U {pz, z) = 
2zpz, or 1 — 2zpz — 2p^ + 2zpl = 0. One easily checks that P{t, z; 1) satisfies the meromorphic 
schema of |141 Thm. IX. 9], and the limit behaviour of X„ + Yn follows. 

The study of the distance between the endpoint and the first diagonal is similar. We first 
refine the enumeration by taking into account the difference X„ — Yn, using a new variable z. 
When estabHshing the functional equation satisfied by T{t,z;u) = T{z;u), one must note that 
walks ending on the right edge of the box have generating function T{t, 1/z; u). The equations 
finally read 

/ tu{l — t'^\\ z t^ 

'i ' ^ ' ^T{z;u) = \~tzT{z;tz) T{z;tz), (7) 



\ (z — tu) [u — tz) I z ~ tu u ^ tz 

P(t,z;u) = T{t, z;u) + T(t, z;u) ^ T{t, z;0) 

with z = 1/z. The kernel method provides a linear equation between T(z; tz), T{z; tz), involving 
the (quadratic) series U{z) = U{t,z) that cancels the kernel of Replacing z by 1/z gives 
a second linear equation, now involving T{z;tz), T{z;tz) and U{z). Thus both series T[z;tz) 
and T[z;tz) can be expressed in terms of U{z) and U{z). It is then straightforward to obtain 
expressions for T{t, z; u) and P[t, z; u) in terms of U{z) and U{z) (we recommend using a formal 
algebra system). In particular, each of these series is algebraic of degree 4. 

By an obvious symmetry argument, E(X„ — F„) = 0. As most continuity theorems involve 
non-negative random variables, we consider n + X„ — Yn, which is clearly non-negative, with 
mean n. The associated generating function is P(tz,z;l). For z in a neighborhood of 1, the 
radius of convergence pz = poi this series is reached when 

p^z* + p5 (p2 + 2 p - 4) z^^ - p3 (4 p2 + 5 p - 5) z-* + p (p - 2 + 5 p2) z2 - 2 p + 1 = 0, 
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and is again a simple pole of P{tz, z;l). One easily checks that we are again in the meromorphic 
schema of [14l Thm. IX. 9], and the limit behaviour of n + Xn — Yn (and consequently, of Xn — Yn) 
follows. B 

4. Enumeration and asymptotic properties of 3-sided prudent walks 
Proposition 7. The generating function of 3-sided walks ending on the top of their box satisfies 



wh 



ere 



1 - tw + t'^ + t^W - J(l - + t-tw + f^w)(l -t-tw - f^w) 

U{w) = U{t;w) = ^ 

is the only power series in t satisfying {U — — tU) = twU{\ — t^), and 



q = q{t) = U [t; 1) = . 

The generating function P{t;u) of 3-sided walks, counted by their length and by the width of the 
box, can be expressed rationally in terms ofU{u) and T{t]u,tu) (see (fTO|) ). When u = 1, this 
gives the length generating function as 

1 / o i^ + t)(2-t-t'^q)\ 1 

P{t;l)^-^—-.ht'qTit;l,t)+^ "'^ 



1 - 2t - 12 I \ ' ' J l-tq / l-t' 

Note that the series U{t] 1) is the algebraic series that occurs in the solution of 2-sided walks 
(Proposition [5]) . 

Proof. In the equation ([3]) satisfied by R{u,w), the only catalytic variable is u (there is no 
occurrence of R{-,w') with w' ^ w). Thus we can apply the standard kernel method: setting 
u ~ U(w) cancels the coefficient of R{u,v), and we are left with 

tu^Rit, H . (^-i^ + tnt^, .)) . (8) 

Recall that T is symmetric in its two (catalytic) variables. In particular, T{tw,w) = T{w,tw). 
The equation ([2]) satisfied by T{u, v) involves the series R{t, u) and R{t, v). We use ([8]) to express 
them in terms of U{u), U{v), T{u,tu) and T{v,tv), and obtain an equation that involves only 
the series T: 

^ _ tuvjl-t^) \ ^ 
(li — tv)(v — tu) J ' 
U(u)~t U(v)-t f tv , ,\ , , f tu , ,\ , s , s 

Now both u and v play catalytic roles. We want to cancel the kernel of this new equation, 
namely the polynomial K{u, v) ~ (u ~ tv){v — tu) — tuv{\ — t^), by an appropriate choice of v. 
Thus V will be a function of u and t. As K{u, v) is homogeneous in u and v, the dependency of 
u in tt is extremely simple: K{u, v) vanishes for v ~ qu, where q = U(t; 1) only depends on t (of 
course, K{u,v) also vanishes for v ~ u/q, but this is not a power series in t). This simplicity is 
crucial to writing the solution in an explicit form. Replacing v by qu in ^ gives 

Observe that tq/{q-t) = jfff = l-\-0{t), while t/(u) = 0{t). Hence ^-U{u) is invertible in 
the ring of power series in t with coefficients in Q[w]. Moreover, j-^^ — U{uq) = 0{t^), so that 
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we can iterate the above equation indefinitely, replacing u by uq, then by uq'^, and so on. The 
net result is the expression of T{u, tu) given in the proposition. 

We now seek an expression for P{t;u), which was given in terms of T{u,u), T(u,0) and 
R{l,u) in the third equation of Lemma [H We wish to express each of these series in terms 
of U{u) (which is known expHcitly) and T{u,tu), which we have just determined. To express 
R{l,u) (or, equivalently, R{l,w)), we combine the case u — 1 of ^ with ([8]). For the other 
two series, namely T{u,u) and T(u,0), we specialize ^ to v — u, and then to v = (using 
U{0) ~ t). Putting together the three pieces gives an expression for P{t; u) in terms of U{u) and 
T{u,tu), which can (for instance) be written as follows: 

2(1 - U{u)){l + t){l - u) / 2^,^, , {t+l)t 



(1 - 2t - t^)ii ^t^tu- t^u) \! + r^muj )-T-f (1°) 

Note that the second term has a factor (u—1): this makes the specialization u = 1 obvious, and 
gives the announced expression of P{t; 1). B 

Where is the group? In the introduction (Section fL2|) . we wrote that a group is associated 
with every Hnear equation with two catalytic variables, and that the group associated with 3- 
sided walks is infinite. What is this group? It is generated by two transformations <& and ^' that 
act on ordered pairs {u,v) and leave K{u,v)/u/v unchanged: 

*(",«)= , *(M,u)=Ut, — 

It is easy to see that they generate an infinite group. Indeed, if one repeatedly applies o $ 
to the pair (m, u), one obtains all pairs [v^"^ lv?'^~^ ^v^"^^^ ju^''') for z > 0. The role of this group 
was first recognized in the study of Markov chains in the quarter plane [13] . See also the more 
recent and combinatorial papers [HI O [6l [HI [TTl [231 [22]. As announced in the introduction, we 
now proceed to show that the series P(^\ 1) is not D-finite. 






Figure 7. Random 3-sided walks of length 400. 



Proposition 8 (Nature of the g.f. and asymptotic properties of 3-sided walks). The 

length generating junction P(t\ 1) of 3-sided walks is meromorphic in the diskV = {t : \t\ < tc}, 
with tc = \/2 — 1. It has infinitely many poles in this disk, and thus cannot be D-finite. On the 
segment (0,tc), the poles are simple and form an increasing sequence {ti)i>o that tends to tc- 

The pole to is the unique dominant singularity of P{t; 1) and coincides with the radius of 
convergence p of the series counting 2-sided walks (Proposition\Bj) . Thus tg = p satisfies l — 2p — 
2p^ + 2p^ = 0, and the number of n-step 3-sided walks is 

Pn ~ K/i" with p — — ~ 2.48 
P 

for some k > 0. 
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Let Wn denote the width of the box of a random n-step 3-sided walk. Then 
where 

1 + p ^ 3p jl + p) (385 - 1148 494 p) 
~ 0.31 and s ~ ^ 5 V ~ 1.41, 



2(1 + 3p) 16(p2 + p-l)(3p- 1)^(1 + 3pr 

and the variable converges in law to a standard normal distribution. 

Proof. We start from the expression of P{t; 1) given in Proposition [71 The polynomial l — 2t — t^ 
does not vanish in V (although it vanishes at tc), and we will show below that q{t) is analytic in 
V, and that 1 — tq{t) does not vanish in this domain. Hence most of our analysis will focus on 
the series T{t; l,t). We will prove that it is meromorphic in V, with a sequence of real positive 
simple poles p ^ ta < ti < t2 < ■ ■ ■ < tc satisfying 

Clearly, our first concern will be the series U{u). 

• The series U{u). The quadratic equation that defines U{u) gives 

U{u)^t + U{u) with Uiu) = -^^^^^, (11) 

which shows that both U(u) and U{u) have coefficients in N[u\. In particular, for all t and u 
such that U{\t\; \u\) converges, U{t\u) also converges and satisfies \U{t;u)\ < U{\t\] \u\). 

For u = 1, we find that the radius of convergence of U{t; 1) = q{t) is at tc. Thus U{t; 1) is 
analytic in V. Moreover, we note that U{tc; 1) = 1. The non-negativity of the coefficients of 
U{t; 1) = q{t) implies that \qit)\ < 1 for t £ X>. 

Consequently, for t G T) and z > 0, the series U{t]q{ty) = U(q^) is (absolutely) convergent, 
and thus analytic in V. Moreover, |tJ7(q*)| < tcU{tc\ 1) ~ tc < 1, so that 1 — tU{q^) does not 
vanish in V. 

• The numerator and denominator of T{t] 1, t). Recall the expression of T(t; it, tu) given in 
Proposition [71 Note that q{ty ^ as i ^ 00, both as a power series in t (because q{t) = 0{t)) 
and for every t £ V (because \q{t)\ < 1). In particular, the term tq/{q — t) — U{q^) = t^ /{q — 
t) — U (g*) converges to t^ /{q — t). Let 



■i>0 



As ^j2^ = 1 + 0{t) and C/(g') = 0{f ), this is a well-defined series in t. We write 



with 
and 

with 



T«;M) = ^ (13) 
JV(() = ^(-l)''iV,(t) 



fc+l 



^) n(T^-W«))n(i-^W))r.w ,14) 



i>k 

U{q'') - t U{q''+^) - t 



t{l-tU{q'')) t{l-tU{q''+^)y 
We will now prove that D{t) and N{t) are analytic in V, so that T(t; l,t) is meromorphic 
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• The series D{t) is analytic in V. As discussed above, every term of the product D{t) is 
analytic in V. We still need to prove that the product converges in V. For t £ V and \u\ < 1, 
the equations ifTTj) imply that \U{t;u)\ < \u\U {\t\; 1) . Consequently, 

\U{q')\<\qitWUm;l), (15) 

with \q{t)\ < 1, so that the series J2i \U{q^)\ is convergent for t G T). The same holds for the 
product D{t). 

• The series N{t) is analytic in T). It follows from the properties of q, U and U that every 
summand Nk (given by lfT4|)) is analytic in P. We will prove the convergence of the series 
J2k>o \-^kit)\ by bounding Nkit). Let us begin with Tk{t). First, we note that 

U{Uu)-t _u{l-t^)U{t;u) 



t{l-tU{t;u)) l-tU{t;u) 

is uniformly bounded (by a constant) for i e I? and |u| < 1: hence Tk{t) is uniformly bounded 
by a constant for t gV. 

Let us now bound the infinite product occurring in Nk{i). We write 



n(i-^^(t;.') 



i>k 



< cxp 

\i>k 



^-^\U{t;q')\ 



The second inequality follows from (fT5|) . 

Hence it suffices to prove the convergence of 



E \M,{t)\ with M,{t) = ( V)'^' n {jh-q ' ^('^'^'0 



Recall that U(t; q{t)^) U{t; 0) = i as fc -> oo. Hence 



Mkit) _q-t { t feA ^ q{q-t) 



Uiq') 



Mk^i{t) <2 yi-tq 'J 1-tq 

Due to the positivity of the coefficients of the series q{t) = U{t; 1) and q{t) — t = U{t; 1), the 
modulus of this ratio is strictly bounded in V by the value it takes at tc, namely 1. Thus 
\Mkit) / Mk-i{t)\ converges, as k grows, to a limit that is less than 1: the convergence of 
J2k \^^k{t)\ follows, and implies that N{t) is analytic in P. 

Consequently, T{t;l,t) = N{t)/D{t) is meromorphic in V. In this disk, ah its singularities 
are poles, found among the zeroes of D{t). The product form (fT2l) of D{t) leads us to study the 
zeroes of each factor, or, equivalently, the zeroes of tq/(q — t) — U(q''). 

• Zeroes of tql{q — t) — U{q^). We fix i > and focus on the real interval (0,<c)- On this 
interval, U{t;q{ty) increases from to 1. The identity 

tq l-tq tU{t; 1) 



q-t 1 - i2 1 _ i2 

shows that tq/{q — t) decreases from 1 to 1/V2 < 1. Thus there exists a unique U in (0,tc) 
such that Uqlu) / [q(ti) - U) = U{q{UY). Moreover, as U{q') > U{q'+^) for t € (0,ic), one has 
tQ < ti < t2 < ■ ■ ■ < tc and 

U{q{Uf)>{) for k>i. (16) 



qiU) - U 
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Each of these zeroes is simple, as the derivative of tq/{q — t) is negative. As T{t] l,t) is mero- 
morphic in V, its poles are isolated, so that the increasing sequence {ti)i can only converge to 

The equation satisfied by p :— to reads q{p) — 2p, which yields the cubic equation of the 
proposition. 

• Existence of infinitely many poles. It remains to prove that each tj, for j > 0, is actually 
a pole of T{t; l,t) = N{t)/D{t), that is, that N{tj) ^ 0. Recall the expression ^ of Nk- For 
t = tj, one has U{q^) = tq/{q — t), or, equivalently, U{q^) — /{q — t). Hence 

N{t) = Y,{-lfNu{t). 
k>j 

We will show the following properties: for t = tj, 

(A) {—ly Nj{t) > 0, and for k > j, the signs of (— l)'^Affe(t) alternate, starting from a positive 
sign: (-l)^-J""iArfc(t) > 0, 

(B) \N,+^it)\ > \N,+2it)\ > \N,+,it)\ > ■ ■ ■ 

The combination of these two properties implies that N{tj) > 0, and in particular that tj is 
indeed a pole of T{t; l,t). 

Let us write t = tj, and study the sign of Nk{t). We begin with the signs of Tk{t) and q — t. 
We have already seen that U{u) — i = U{u) > and tU{u) < 1 for u G (0, 1). In particular, 

Tk{t)>l and q-t>0. 

Let us move to the infinite product occurring in Nk{t). Now for i > k> j, 

U{q=) 

as g e (0, 1) and U{u) is an increasing function of u, for < u < 1. 

We are left with the sign of t/ (1 — tq) — U (g*^^), for i > 0. We will prove that, still denoting 

t = tj, 

U{q'+') < < Uiq^+'). (17) 

Given that the sequence J7(g') decreases as i grows, this gives 

— U(q'+^) <0 forO<i<j, 

1 — tq 

> for i > j + 1, (18) 

and concludes the proof of Property (A). 

The key to prove lfT7|) is to observe that the function u^U{t; u) is increasing for u G (0, 1), 
with an explicit inverse, easily derived from the quadratic equation defining U{u): 

(U{u)-t)(l-tU(u)) 
tU{u){l-f^) 

Thus (fTT]) is equivalent to 



J+2 . q{l-tq-t') 



Recall that we write t = tj. Hence U{q^) ~ tq/ [q — t), so that 

^ _ jq-t- t^q) 



q{l-t^){q-ty 



Hence we are left with proving that 



q{q-t- t^q) ^ q{l - tq ~ t') ^ {q - t - t'q) 
q — t 1 — tq Q — t ' 



PRUDENT SELF-AVOIDING WALKS 



17 



which is easily seen to hold: the first inequality boils down to q < 1, and the second one to 
q — t — t^q — t^q > 0. But the quadratic equation defining q gives 



q - t - t^q - t-^q = tq{l + t){l - 
It remains to prove Property (B), that is, for k > j + 2, 



2t) > 0. 



\Nkit)\ 



\Tkit)\ 



^ - uiq'^: 



< 1. 



First, note that the sequence [U{q''))k is decreasing, 
it suffices to show that for fc > j + 2, 



t 



Uiq") 



l~tq 

Given lfT6|) and ifTSj) . we are left with proving 

t 



< 



tq 



\Tk-i{t)\ 

Hence the same holds for [Tk[t))k. Thus 
- U{q'^) 



1 - tq 



< 



tq 



t' 



which boils down to q^ < 1 and is immediate. 

This concludes the proof of the properties of T{t; l,t) stated at the beginning of the proof. 
These properties also hold for P(t; 1), as shown by the last equation of Proposition [71 

• The asymptotic number of 3-sided walks. We have proved that T{t; 1, t) and P{t; 1) are 
meromorphic in V, with a smallest real positive pole at to = p. By Pringsheim's theorem, p 
is the radius of convergence of P{t; 1). In order to prove that the coefficients of P{t; 1) behave 
asymptotically as Kp~^ , we need to prove that there exists e > such that P{t] 1) has no 
singularity other than p in the disk of radius p + e. 

As P(i; 1) is meromorphic in V, its poles are isolated. Thus we only have to prove that there 
is no other pole of modulus p. Recall that all poles t of P(t; 1) are roots of tq/{q — t) = U{q^) 
for some z > 0. Let t be a complex number of modulus p, with t ^ p. The positivity of the 
coefficients of Uit] u) imply that for i > 0, 



\U{t-q^)\<U{p-q{pr)<U{p;l) 



while 



tq 



1 - 



tij{t]i) 



> 1 - 



pU{p- 1) 



U{p;l) 



Hence < \tq/{q — t)\ for all i and t cannot be a pole. 



• The width of the box. The series P{t\ u) that counts 3-sided prudent walks by their length 
and the width of the box is given by ifTO]) . A perturbation of the singularity analysis that we 
have performed for P{t\ 1) reveals that, for u in some neighborhood of 1, P(t; u) admits a unique 
dominant singularity, which is a simple isolated pole obtained when U{u) = tq/(q — t), or 



t {t - If (t + iyu'^ + {t^ + 2t'^-2t+ 1) (t 



1)' it 



t{l-2t^ -t + t^ + 2t^ 



We are again in the meromorphic schema of [HI Thm. IX. 9], and the limit behaviour of Wn 
follows. ■ 



5. Enumeration and asymptotic properties of triangular prudent 

WALKS 

We now turn our attention to the triangular prudent walks of Fig. [3l Recall that the box of 
such a walk is a triangle that points North. If this box has size k, we say that the walk spans a 
box of size k. 
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Proposition 9. The number of triangular prudent walks that span a box of size k is 

Pfc = 2^-i(fc + l)(fc + 2)!. 

More precisely, the number of triangular prudent walks that span a box of size k, and end on the 
right edge of their box at distance i from the North corner ( and thus at distance j ^ k — i from 
the South- West corner) is 

2Hfc + 2)! „ . „ 
«/ « = or .7 = 0, 

2''{k + 2)\ 



6 



otherwise. 



Proof. We specialize the functional equations of Lemma [4] to t = 1 . Remarkably, the kernel 
of (HJ reduces to 1. Denoting R{u,v) = R{l-u,v), we have 

v) = l + 2u u) + 2v R{v, v). 

V ^ U U — V 

This deprives us of our favourite tool (how could we cancel a kernel that reduces to 1?), but still 
yields a simple solution. Let Rk{u,v) be the homogeneous component of degree k in R(u^v). 
With the notation fij used in the proposition, 

Rk{u,v) = ^ h^ju'-v^. 

i+j—k 

The functional equation satisfied by R{u,v) is equivalent to the following recurrence of order 1: 

•- u{v -2u)Rk-i{u,u) - v{u-2v)Rk-i{v,v) 

Rk[u,v) = 2 , 

V — u 

with initial condition Ro{u,v) = 1. One readily checks that the polynomial 

2'=(fc + 2)! / f, „fe+i_„fc+i 



Rk{u,v) = — [u^+v 



U ' — V 



satisfies this recursion, and yields the values fij given in the proposition. The expression of pk 
then follows using the second equation of Lemma HI ■ 

Remark. The fact that the numbers pk grow faster than exponentially is not unexpected. In 
particular, given a square of size k, the number of partially directed (1-sided) walks of height k 
that fit in this square is (fc + 1)'^+^. Indeed, such walks are completely determined by choosing 
the abscissas of the k vertical steps, of the starting point and of the endpoint. 
We now move to the length enumeration of triangular prudent walks. 

Proposition 10. Set u = {i+tx)(i+t'^x) ' '^'^^re x is a new variable. The generating function of 
triangular prudent walks ending on the right edge of their box satisfies 

. , ^. 2^ V-*^'^'U^*(1-2<'))V Xt^ \ 

Rit; uM .^t)(i + xt^) g [i^-^Vk ^''^ 

where we have used the standard notation 

{a;q)„ = (1 - a) (1 - ag) • • • (1 - ag"-i). 

The generating function of triangular prudent walks, counted by their length and the size of the 
box, is 

") = 1 + i_t-2fa(l + t) (1 + ^ (2«(1 + - 1) i?(t; u, tu)). 
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Comments 

1. The parametrisation of u in terms of the length variable t and another variable x is just a 
convenient way to write R{t]u,tu). Equivalently, we have 



R(t;u,tu) = (l + Xt){l + Xt^) 



where 



X = X{u) = 



\-t-ut-ut^ - - t){l -t-2ut- 2ut^ + uH"^ - uH^) 

2M^ 



is the only power series in t satisfying u(l + tX){\ + t^X) = X{1 — t). In particular, the length 
generating function of triangular prudent walks is 

&t{l + t) 



P{t-l) 



I -it- 2t^ 



{l + t{l + 2t)R{t-l,t)) 



where 



and 



Note that 



{Y{l-2t^)f ( 



-2t2' 



Y = tX[\) 



\-2t-e- ^(l-i)(l-3t-t2-t3) 



Y 



1 - t 



2*2 



(i + y)(i + tr). 



2. The value of R{u,tu) is especially simple when x = 1/(1 — 2*2)^ that is, for u 
Indeed, for this value of u, 

, (l + t-2t2)(l-<2) .fc+C^r)/ <3 X 



(20) 
(21) 

(22) 
(23) 

l — 2t^ 
(l-t^)(l+2t) • 



^(-n.-n(:-^)) 



The product form follows from the following identity [H Corollary 2.7]: 

m>l 



E 

n>0 



{t\t)r 



(24) 



This can be used to prove that neither R{t; u, tu) nor P(i; u) are D-finite. However, we will 
derive from l(20H2T|) the following finer result on the length generating function P{t; 1). 



> 



Figure 8. Random triangular prudent walks of length 500. 
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Proposition 11 (Nature of the g.f. and asymptotic properties of triangular prudent 
walks). The length generating function P{t] 1) of triangular prudent walks is meromorphic in 
the domain V = {t : \t\ < 1} \ [tc, 1), where tc ~ 0.295... is the real root of 1 - M - t^ - . In 
this domain, it has infinitely many poles, so that it cannot he D-finite. These poles accumulate 
on a portion of the unit circle, so that P{t; 1) has a natural boundary. 

The series P{t; 1) has a unique dominant singularity, which is a simple pole p = (\/l7— 3)/4 ~ 
0.280. Hence the number of triangular prudent walks of length n satisfies 

Pri ~ K I — ^ 

for some positive constant k. 

Let Sn denote the size of the box of a random n-step prudent walk. Then the mean and 
variance of Sn satisfy 

where 

\ 1 \ .2 12 



17 



m = - 1 H ■= and s 



2 V \/T7/ 17^/T7 
and the variable '^^y^" converges in law to a standard normal distribution. 

Proof of Proposition [TOl The kernel of the functional equation ^ reads: 

K{u, v) = {u — tv){v — tu) — tuv{l — + v). 

Again we want to cancel it by an appropriate choice of u and v. The kernel is not homogeneous 
in u and v, as it was in the case of 3-sided walks, but it has another interesting property: the 
curve K{u, v) = has genus 0, and thus admits a rational parametrisation, namely 

KiUi.)Mt^))^0 for ^W- (,4;~;^2,) . 

Now take an indeterminate x, and set u ~ U{x) and v = U{tx). The series R{u,v) is well- 
defined (as a series in t with coefficients in R[a::]). Since the kernel vanishes, it follows from |[4]) 
that 

, V — 2tu „ , , , s It — 2tv „ , 

1 + tu(l + t) R{u, tu) + tvil + t) R{tv, v) = 0. 

V — tu u — tv 

Recall that R{u, v) ~ R{v, u) for all u and v. Denoting $(a;) = R{u{x), tu{x)), this equation can 

be rewritten as 

{1 + xt){l + xt"^) xt'^{l+xt){l-2t'^ - xt^) 
^""^ " l-xt(l-2<2) + {l + xt3){l-xt{l~2t^)y^''^'' 

Iterating it gives the value (fT9l) of R{u,tu). 

The expression of P{t;u) in terms of R{u,tu) follows from |(5|), after specializing ^ to v = u 
and then v ~ 0. B 

Where is the group? We have solved a new linear equation with two catalytic variables u 
and V. Again, two transformations $ and ^' leave K{u,v)/u/v unchanged: 

„,2 



$(u,z;)= — — -pr-^,v], *(m,u)= u, 



^u{l + v — t^v)' / ' \ ' 7j(l + u — i^u) ^ 

They generate an infinite group. Indeed, if one repeatedly applies 'I'off) to the pair {U{x), U{tx)), 
one obtains all pairs {U{t^^x),U{JP'^'^^x)) for z > 0. We now proceed to show that the series 
P{t] 1) is not D-finite. 

Proof of Proposition [TTl Recall the expression ((20|) of P{t] 1). The denominator 1 — 3i — 2t^ 
vanishes at p and will be responsible for the simple dominant pole of P{t; 1). We will see that 
the series R{t;l,t) = R{l,t), given by (f2T|) . has a larger radius of convergence than P(t;l). 
More precisely, we will prove that R{l,t) is meromorphic in T>, with a unique real pole at 
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to — 0.288 e {p,tc), where 4to'* + 2to'^ — 6to^ — 2<o + 1 — 0, and infinitely many non-real poles 
that accumulate on a portion of the unit circle. 

• The series Y{t). The square root occurring in the series Y given by l(22|) vanishes at t = 1, 
at t = tc — 0.295 and at two other points of modulus 1.83... > 1. Hence Y has radius tc, but 
can be analytically continued in the domain V. In this domain, Y{t) is bounded. Moreover, it 
is easily seen to be increasing on the segment (— l,tc)- 

• The series is meromorphic in V. Let us write 

where 
and 

YP 



N{t) = (1 + + tY) J2 {Y{1 - 2t^))' (^^'^'(1 - 2^'); 

fe>o ^ ^ 

Clearly -/V(t) and D{t) are analytic in I? (as \t\ < 1). The cancellation of (1 — 2t^) does not 
create any singularity, as 

(1-2^2)'^ (^-K^-t^ ^ [\-2e -Ye)[\-2t^ ~Yt^)---{\-2e ~Yt^"-^). 

Hence is meromorphic in V, and its poles are found among the values of t such that 

Y(i:)t\\ ~ 2t^) = 1 for some £ £ N. All poles are simple, as Y{t)t\l - 2F) = 1 implies 
y(t)t™(l - 2t2) = t™-^ ^ 1 for m 7^ 

• Real poles of R{l,t). A standard study of the functions t ^ Y{t) and t ^ l/t7(l ^ 2t2) 
on the interval (— l,tc) reveals that the only possible real pole of is at to — 0.288 < tc, 
where y(to) = 1/(1 - 2tg). This implies 4 V + 2 to^ ~ 6 to^ - 2 <o + 1 = 0. 

• Existence of infinitely many poles. Let us call a critical value any t £ V satisfying 
Y{t)t^{l — 2t^) = 1 for some £. We have just exhibited a real critical value t^. Fig. [9] shows 
the critical values obtained for various values of £. We will first prove that V contains infinitely 
many critical values. Then, we will show that almost of them, and in particular to, are poles of 
R{l,t). 

To prove the former point, we prove that there exist < 9o < 9i < n such that every point 
e'^ with Oq < \6\ < Oi is an accumulation point of critical values. As Y{z) = Y{z), it suffices to 
consider the case 9 > Q. Let 9 S (0, tt) be of the form pn/q for two integers p and q. Let I = 2kq 
with fc > 0. Take t in the vicinity of e**, that is t = e**(l + s) with s small. As Y is analytic in 
the neighbourhood of e**, 

Y{t)t\l - 2t') = r(e'^)(l - 26^'") exp(£s + 0{es^)){l + 0{s)). 




Figure 9. The roots of r(t)t7l-2t2) = 1 lying inside the domain I? for ^ 20, 
i = 30,£ = 200. 
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This shows that there exists, in the neig hbourhood of e*^ a solution t of Y{t)t'^{l - 2t^) = 1 for 
every (large) £ = 2kq. This solution reads t = e**(l + s) with 

s~-ilog((l-2e2''^)r(e'^)). 

However, t will only be critical if it lies in V, that is, if |t| < 1. But 

|ip==l-^log(|l-2e2''^||y(e*«)|)+o(ri), 

so that t is critical (for £ large) if and only if f{e) |1 - 2e'^''^\\Y{e''^)\ > 1. The function / is 
continuous on (0,7r). Moreover, at t = e'''/^^ one has (1 - 2t'^) = 3 and Y{t) = -l + i± e~"/'', 
so that /(7r/2) > 1. This proves the existence of an interval (6*0, 611) where every point 6 satisfies 
f{6) > 1 and is thus an accumulation point of critical values. Fig. [TO] shows a plot of f{0), and 
gives estimates for the best possible values 9o ~ 0.59 and di ~ 2.28, in good agreement with 
Fig.H 



Figure 10. The function f{e) |1 - 2e^'''\\Y{e''')\. 
LetteV satisfy Y{t)f{l - 2t^) = 1. Then t is a pole of R{l,t) if and only if N{t) ^ 0. But 

N{t) = (1 + r)(i + tY) tr.it-^^ { (SL^^ 'A (^''■'"'; 

fc>0 I Jk 

If fc < £, the fcth summand is zero because of the term {t^^'^~^\i)oa- We thus take k ~ 1 + j, 
with j > 0. We also rewrite (1 + Y){1 + tY) using This gives 



j>0 

,-i-(T)2^f .t\ (ft) n(i+t™)A-^^ 

l-2t^\{l~2t^)^')^^'>°°^^^^ >\ (l-2t2)2 



m>l 

The final product expression relies on l(24|) . 

Could this product vanish for our critical value t G 2?? As |t| < 1, the only factors that 
might be zero are of the form 1 — P /{I — 2t^)^, for j >2 — £. Each such factor has only a finite 
number of zeroes in V. In particular, those for which j < 0, taken together, will only vanish 
for a finite number of critical values, while there are infinitely many such values. For j > 0, 
the roots of 1 — t-' /(I — 2^^)^ that lie in the unit disk satisfy 11 — 2^^1 < 1. That is, they are 
inside the 'glasses' curve |1 — 2t'^\ = 1 plotted on Fig.[TTJ In particular, they remain away from 
the accumulation points we have exhibited on the unit circle, so that only finitely many critical 
values are non-poles. 

Could to, the unique critical value found on the segment (— l,tc), not be a pole? For t — to, 
one has 

(1-2*2)2 > (l_2t2)2 > 1 for j < 0: 

(Tzft2)2 < ^T=yr - 0-4 < 1 for J > 1- 
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Figure 11. The singular landscape oi P{t; 1): this series is meromorphic inside 
the unit disk except on a cut [tc, 1), with tc ^ 0.295. There are two simple poles 
on the real axis before tc, namely p c=i 0.280 and to c=i 0.288. Almost all critical 
values, which accumulate on a portion of the unit circle, are simple poles. The 
'glasses' curve is |1 - 2i2| = 1. To the right, a zoom in the neiffhbourhood of 



Hence N{to) ^ 0, so that to is indeed a pole of R{1, t), and, by Pringsheim's theorem, its radius 
of convergence. 

• The singular landscape of R{l,t) and P{t; 1). We can now conclude as to the nature of 
the series i?(l,i): this series is meromorphic in the domain V, with a dominant pole at and 
infinitely many poles that accumulate on a portion of the unit circle, {e'^ : Oq < \9\ < 6i}. 

In view of (f20| . the same holds for the series P{t; 1), with to replaced by p as the denominator 
1 — 3t — 2t^ induces a new pole p < to. This is illustrated in Fig. [TTJ Clearly, P{t; 1) has no 
other pole of modulus less that to- Hence the number of n-step triangular prudent walks is, as 
announced, asymptotic to np^". 

• The size of the box. Consider the expression of P{t; u) given in Proposition [Tol The 
singular pattern we have found for u = 1 persists when u is in a neighborhood of 1. The radius 
of R{u,tu) is reached when Xt{l — 2t^) =■ 1, but the pole of P(t;u) given by the cancellation 
of 1 — t — 2tu{l + t) is smaller in modulus. We are again in the meromorphic schema of [14t 
Thm. IX. 9], and the limit behaviour of Sn follows. 



6.1. Generating trees and random generation 

A generating tree is a rooted tree with labelled nodes satisfying the following property: if two 
nodes have the same label, the multisets of labels of their children are the same. Given a label 
£, the multiset of labels of the children of a node labelled £ is denoted Ch(£). The rule that gives 
Ch{£) as a function of £ is called the rewriting rule of the tree. 

Consider a class C of walks closed by taking prefixes: if w belongs to C, then so do all prefixes 
of w (this is the case for our classes of prudent walks) . One can then display the elements of C 
as the nodes of a generating tree: the empty walk sits at the root of the tree, and the children 
of the node labelled by w are the walks of the form ws belonging to C, where s is a step. Fig. [12] 
(top) shows the first few levels of the generating tree of 2-sided walks. An even simpler example 
is the tree of directed walks (walks with N and E steps), which has the following rewriting rule: 



In this context, the (uniform) random generation of a walk of C of size n is equivalent to the 
random generation of a path of length n in the tree, starting from the root. This path can be 
built recursively as follows: one starts from the empty walk (the root of the tree), and builds 



0.28. 



6. Random generation 
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Figure 12. The generating tree of 2-sided walks (top) and one of its good 
labellings (bottom). 



Ex(u', m) = 



the walk step by step, choosing each new step with the right probability. If at some point one 
has obtained a walk w of length n — m, with to > 0, then the probability of choosing s as the 
next step must be 

^, , , Ex(it;s, TO — 1) 

F(s\w) = ^ / — ^ 

Ex(i/7, to) 

where Ex(?i',TO) is the number of extensions of w of length to, that is, of walks w' of length m 
such that the concatenation ww' belongs to C. Clearly, Ex(w,m) = if w ^ C. For w G C, the 
numbers Ex(u', to) can be computed inductively: 

1 if m = 0, 

T,wsech{w) Ex(ws, m - 1) otherwise. 

This is the basic principle of recursive random generation, initiated in [21]. The (time and space) 
complexity of the procedure that generates an 7i-step walk of C depends on how efHciently one 
solves the following two problems: 

(1) compute the numbers Ex{w, to), for all walks w of length n — m, 

(2) given a walk w of length less than n, decide which steps s are admissible, that is, are 
such that ws belongs to Ch(?u). 

We will explain how both tasks can be achieved by introducing good labellings of the class C. A 
function L defined on C is a good labeUing if for all w & C, the multiset {L{w') : w' S Ch{w)} 
depends only on L{w). This means that the tree obtained by replacing each node w of the 
generating tree by L{w) is itself a generating tree. We will illustrate this discussion with the 
class of 1-sided walks. 

(1) The numbers Ex(?i',TO). When C is the class of 1-sided walks, these numbers only depend 
on the direction (horizontal or vertical, H or V for short) of the final step of w (and of course 
on m). Denote this direction by L(w). Then it is easy to see that L is a good labelling: the 
tree obtained by replacing each node w of the generating tree by L{w) is itself a generating tree, 
with rewriting rule 

(By convention, the label of the empty walk is set to V.) 

For a general class C, any good labelling L allows us to compute the numbers Ex(w,to). 
Indeed, they can be rewritten as Ex(L(ui),to), with 

fl ifTO = 0, 

Ex(£, to) = I J2e'ecm Ex(f , m - 1) otherwise. ^^^^ 
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In the case of partially directed walks, we have 

Ex{V,m) = Ex{V,m- 1) + 2Ex(i?,m - 1), Ex{H,m) = Ex{V,m- l)+Ex{H,m - 1). 

A good labelling is compact if it takes its values in Z'' , with k fixed. In this case, it is usually 
easy to determine Ch{£), given £. In order to fasten the computation of the extension numbers, 
one tries to minimize the number of distinct labels L{w) occurring in the first n levels of the 
generating tree. This solves Problem (1) above. 

(2) Admissible steps. We are left with the second problem: determine the admissible steps. 
Moreover, we also need to decide, for every admissible step s, which of the elements of Ch(L(u;)) 
is the label of ws. That is, we need a way to correlate new steps with new values of the label. 
In our random generation procedures, this will be achieved thanks to a second good labelling P 
having the following properties: 

- it refines the first one; that is, P{w) = {L{w)\M{w)), and the rewriting rule of P extends 
that of L, 

- the last step of w is S{P{w)), for a simple function S. 

For instance, in the case of 1-sided walks, we keep the label V for walks ending with a vertical 
step, but introduce two versions of H, namely (H \ 1) (for a final E step) and {H | — 1) (for a 
final W step), and refine the rule ((25| as follows: 



The last step function is given by S{V) = N, S{H \ 1) = E, S{H = W. This new 

rewriting rule is less compact than ((25l) . and not needed for the calculation of the extension 
numbers (why compute and store separately Ex((i/ | 1), to) and Ex((i/ | — 1), to) while they are 
equal?). However, it is needed for the generation stage. 

In our random generation procedures below we use the i-labels to compute the extension 
numbers, and the P-labels to determine the admissible steps and associate with them the correct 
new label. Accordingly, the algorithm involves two procedures: 

• The procedure Ex{£,m) computes the extension numbers using i(26|) . 

• The generation procedure reads (denoting p[l] = £ for a P-label p — {i\mj): 

gen:=proc(n) 

Initialisation w: = [] : p:=P(w): f:=p[l]: m:=n: 
while m>0 do 



[pi, - • ■, Pfc] := Ch(p) : 

for i from 1 to k do ii:='Pi[11 od: 

U:=Uniform(0,l) : d:=U*Ex(£,m) : 

i:=l: 

while d> Ex(£i,m-1) +■ ■ • +Ex(£j ,m-l) do i:=i+l: od: 
p:=Pi: s:=S(p): w:=ws: 



Since all our procedures have this form, the following sections only describe the labellings L 
and P that we use for each class of walks. In all cases, the number of children of a P-label is 
bounded, so that the generation procedure is achieved in linear time. We first discuss the simple 
case of 2-sided walks, and then the more complex case of general prudent walks (on the square 
lattice). Once the principles underlying the corresponding labellings are understood, it is not 
difficult to adapt them to 3-sided prudent walks and to triangular prudent walks. 

6.2. Two-sided prudent walks 

A non-empty 2-sided walk w can be of three different types, /, C or F, depending on the 
nature of its final step s: 

(/) either s is an inflating step, which moves the North or East edge of the box; 




od: 
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(C) or s walks along the North or East edge of the box, moving the endpoint closer to the 
NE corner, 

(F) or s walks along the North or East edge, moving the endpoint further away from the NE 
corner. 

We label w by L{w) = (T; i), where T is the type of w, and i is the distance of the endpoint of 
w to the NE corner of its box. It is not hard to realize that L is a good labelling: the labels of 
the children of a node w of the generating tree only depend on L{w). These labels are described 
by the following rewriting rule: 

[ (/;0) if 1^0 ^ ^ '"^ " 

We do not define the label of the empty walk. Two of its four children have label (/; 0) (cor- 
responding to North and East steps), while the other two (corresponding to South and West 
steps) have label {F; 1). The top of the label based tree is shown in Fig. [12] (bottom) . 

The number of extensions of length m of a walk labelled £ can now be computed using l(26|) . 
For a walk of length m and label {T;i), one has i < m. Hence 0{n) different labels occur in 
the first n levels of the tree, and the calculation of the extension numbers Ex(£, to) needed to 
generate a random walk of length n requires O(n^) operations. 

The L-label L{w) does not determine the last step s of w. For the generation stage, we 
refine it by defining P(w) = {L{w) \ s). The last step function is of course S{P{w)) = s. The 
generating tree based on the labelling P has rewriting rule: 

^) 

3-s) 

1 - s) if i > ' 
1 - s) if i 

(/;z I 3-s) 

{F;i + 1 I s) ■ 




{C;i\ s)^ { iC;i-l \ s) if i > , {F;i 



We have written s = (resp. 1, 2, 3) for a North (resp. East, South, West) step. The four 
children of the root have labels (/; | 0), (J; | 1), (F; 1 | 2), (F; 1 | 3). 

Remark. As shown in |llj . the language on the alphabet {N, S, E, W} that naturally encodes 2- 
sided walks is non- ambiguous context-free. This implies that one can generate these walks using 
another type of recursive approach, based on factorisations of the walks [TSj. This approach 
achieves a better complexity than our naive step-by-step construction: only 0(n) operations are 
needed in the precalculation stage. Moreover, it can be optimized further using the principles 
of Boltzmann approximate- size sampling |12j . As the other classes of prudent walks we have 
studied cannot be described by context-free grammars, we have focussed above on the step-by- 
step recursive approach. Independently, we have also implemented the Boltzmann approach, 
which produced the second and third walks of Fig. [6l 

6.3. General prudent walks on the square lattice 

A non-empty prudent walk w can be of two different types, / or A, depending on the nature 
of its final step s: 

(/) either s is an inflating step, which moves one of the edges of the box, 

[A) or s walks along one of the edges of the box. 
Consider the last edge of the box that has moved while constructing w. Say it is the North 
edge (the other cases are treated similarly after a rotation). Then w ends on this edge. If the 
final step s is infiating (North), or walks along the top edge in clockwise direction (East), let i 
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(resp. j) be the distance between the endpoint of w and the NE (resp. NW) corner of the box. 
Otherwise, s walks along the top edge in a counterclockwise direction (West): exchange the roles 
of i and j. The top edge has length i + j. Let h the other dimension of the box. We label w by 
L{w) ~ (/; {«, j}, h) if w is inflating, by L{w) = {A; h) otherwise. If t = j , the set has 
to be understood as the multiset where i is repeated twice. 

It is not hard to reahze that the labels of the children of a node w only depend on L{'w). 
These labels are described by the left part of the rewriting rules given below (ignoring for the 
moment what follows the vertical bar). As i,i and h are bounded by n, the calculation of the 
extension numbers needed to generate a random walk of length n will require O(n^) operations. 

The label L{w) does not determine the last step s of w. For the generating stage, we refine 
it into a P-label as follows: 

- for a walk of type /, we keep track of the ordered pair (i, j) and of the last step s; hence 
the P-label reads (/; {i, j}, h \ s), 

- for a walk of type A, we keep track of the last step s and of its direction d around the 
box (d = 1 if s goes in a clockwise direction around the box, d = —1 otherwise). The 
refined label thus reads (A; i,j,h\s,d). 

The combined rewriting rule reads: 



(-^;{«,j},/i I «,i,s) 



(/;{^,J■},/^+l 
{A;i - l,j + I, h 
{A;j -l,i+l,h 
(/;{0,/i},j + l 
{I;{0,h},i+1 



hJ,s) 
s + 1, 1) 
s- 1,-1) 
/i,0,s + l) 
0,/i, s- 1) 



if i > 
if j > 
if i = 
if J-0 



with 



d) 



{I-{t,j},h+l 

{A;i,j,h\s,d)^{ {A;i - l,j + l,h \ s,d) if i > , 

(/;{0,/i},i + 1 I h-^0,0-^h,s) if i = 



d 



(27) 



i if d = 1 , 
j ifd=-l. 

Again, we have written s = (resp. 1, 2, 3) if s goes North (resp. East, South, West), and 
these values are taken modulo 4. The four children of the root have labels (/; {0,0}, 1 | 0,0, s) 
for s = 0,1,2,3. 



6.4. Three-sided prudent walks 

A non-empty 3-sided walk can be of five different types, /„, Ih, A, C, or F, depending on the 
nature of its final step s: 

(J„) either s is a vertical infiating step (it moves the North edge of the box), 

(Ih) or s is a horizontal infiating step (it moves the East or West edge of the box), 

{A) or s walks along the top edge, 

(C) or s walks upwards on a vertical edge, moving the endpoint of the walk closer to the top 
end of the edge, 

(F) or s walks downwards on a vertical edge, moving the endpoint further away from the 
top end of the edge. 

For a walk w of type /„ or A, we denote by i and j the distances between the endpoint of 
the walk and the NE and NW corners, with the same convention as for general prudent walks 
(Section l6.3p . Otherwise, i denotes the distance between the endpoint of w and the top end 
of the vertical edge where w ends, and j denotes the width of the box. The label of a walk w 
of type T is L{w) — if T = and {T;i, j) otherwise. Then L is a good labelling, 

described by the first part of the rewriting rules given below. The calculation of the extension 
numbers needed to generate a random walk of length n requires 0{n^) operations. 
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The L-labelling does not determine the last step s of w. For the generating stage, we refine 
it as follows: 

- for a walk of type /„, we keep track of the ordered pair (i, j), 

- for a walk of type A or Ih, we keep track of the last step s, 

- for a walk of type C or F, we keep track of the direction d of the last step around the 
box: d = 1 if the last step goes in a clockwise direction around the box or consists only 
of South steps, d = —1 otherwise. 

The last step is explicitly given by the refined label for types A and Ih , and is otherwise obtained 
via the last step function S defined by S'(/„; ■ ■ •) = S{C; . . .) = 0, S{F; . . .) — 2. The generating 
rules read, with the notation l(27|) : 



{Iv;{i,j} I i,j 



{Iv;{i,j} 
{A;i + l 
{A;j -1,1+1 
ih;0,j + l 
ih\0,i + l 



1) 
3) 
1) 
3) 



if i > 
if j > 

if j = 



{lh;i,j I s) 



ih;ij + 1 
{F;i + l,j 
{C;t-l,j 



2-s) 

s-2) if i > 



j— 0,0^j) if; 







(A; I, J I s) 



{Iv;{i,j} 
{A;i-l,j + l 



s) ifi>0 
s) if i = 



r ilh;ij + l I 2 + d) 

(C;i,i \d) ^ I iC;i-l,j I d) ifi>0 

I (^;{0,j} I J-0,0^j) ifi = 



r {F;i + l,j I d) 
{F;i,j \d) ^ { ih;i,J + l I 2-d) 

[ ih■,^^ I 3) ifj^O 



The four children of the root have labels (/,,;0,1|1), (/,j;0,l|3), (/„; {0, 0}|0, 0) and (i^;l,0|l). 



6.5. Triangular prudent walks 

The principles that underHe the generation of triangular prudent walks are the same as for 
prudent walks on the square lattice (Section 16. 3p . The type of a walk w depends on whether 
its last step s is inflating (type /) or walks along the box (type A). Consider the last edge of 
the box that has changed while constructing w. Assume it is the right edge (the other cases 
are treated similarly after a rotation). If s points up (NW or NE), then i (resp. j) denotes the 
distance between the endpoint of w and the top vertex (resp. bottom right vertex) of the box. 
Otherwise, the roles of i and j are exchanged. Define L{w) = {T; where T is the type of w. 
Then L is a good labeUing. From this, one computes the extension numbers needed to generate 
a walk of length n in time 0(n^). 

During the generation stage, we also keep track of the last step s of w, and of its direction 
d around the box: if the last edge that has changed is the right one, then c? = 1 if s goes East 
of South-East, d = —1 otherwise. The other cases are treated similarly after a rotation. The 



PRUDENT SELF-AVOIDING WALKS 



29 



combined rewriting rules read: 



I s,d) 



{A]i,j I s,d) 



Again, the steps are encoded by 0, 1, . 
NW step is encoded by 0. The six children of the root have labels (/; 0, 1 
(0,-l),(l,l),(2,-l),(3,l),(4,-l),(5,l). 





+ 1 




s,d) 






1 1 

r 1 




s — a, —a ) 




I A - n 


1 T 

1, z - 


- i 


s — ZcL, — a j 




( A. ;^ 

\A, I - 


1 


- i 


s + a,a) 


11 I > V 




-1 

+ i 




s + a, — uj 


It I = 


II- i 1 






s + 2d d) 


Hi = 




+ 1 




s — d,d) 




{I; i 


+ 1 




s - 2d, -d) 




{A-i- 


1,3- 




s,d) 


if i > . 




+ 1 




s,-d) 


if i = 








s + d,d) 


if i = 


...,5 


m clockwise order, and considered modulo 6. The 



s, d), with (s, d) 



7. Some questions and perspectives 

The enumeration of general prudent walks remains an open problem, although we have es- 
tablished a functional equation defining their generating function (Lemma [3]). Recall that this 
g.f. is expected to be non-D-finite, and to have the same radius of convergence as the g.f.s of 
2-sided and 3-sided walks (Propositions [6] and [HI) [TO] , 

It may be worth investigating the number of prudent walks that span a given box. This 
number is indeed simple in the case of triangular prudent walks (Proposition [9]) , even though 
the length generating function is not D-finite. 

Regarding random generation, we think that we have optimized the step-by-step recursive 
approach. However, it only provides walks with a few hundred steps. We believe that a better 
understanding of our results, and of the combinatorics of prudent walks, should lead to more 
efficient sampHng procedures. 

All the walks we have been able to count move away from the origin at a positive speed. 
Would it be possible to introduce a non-uniform distribution on these walks that would make 
them significantly more compact? Recall that the end-to-end distance in an n-step self-avoiding 
walk is expected to grow Hke n^^^ only. 

In the same spirit, it should be possible to study kinetic versions of prudent walks, where at 
each tick of the clock, the walk takes, with equal probability, one of the available prudent steps. 
The resulting distribution on n-step walks is no longer uniform. This version is actually the one 
that was studied in [30] and [27] . Clearly, the random generation of these kinetic walks can be 
performed in linear time, and does not require any precomputation (in terms of generating trees, 
one chooses uniformly one of the children). Random trajectories differ significantly from the 
uniform case (Fig. [13]). Their behaviour has very recently been explained by Beffara et al. [4]. 




Figure 13. Kinetic prudent walks with 1000, 10000 and 20000 steps. 
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